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ABSTRACT 

Receiver timing synchronization of an optical Pulse-Position 
Modulation (PPM) communication system can be achieved using a phase-locked 
loop (PLL), provided the photodetector output is suitably processed. The 
magnitude of the PLL phase error is a good indicator of the timing error at 
the receiver decoder. The statistics of the phase error are investigated 
while varying several key system parameters such as PPM order, signal and 
background strengths and PLL bandwidth. A practical optical communication 
system utilizing a laser diode transmitter and an avalanche photodiode in 
the receiver is described, and the sampled phase error data are presented. 

A linear regression analysis is applied to the data to obtain estimates of 
the relational constants involving the phase error variance and incident 
signal power. 
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CHAPTER 1. INTRODUCTION 

Several techniques have been developed to establish effective inter- 
satellite communication systems. Among the more promising of these is the 
optical communications link. The design and characteristics of optical 
Pulse-Position Modulation (PPM) communication systems have been treated by 
several researchers [1] - [2]. These studies have demonstrated that the PPM 
signaling format offers several advantages over other digital signaling 
schemes when employed in free space optical links. 

Pulse-Position Modulation refers to the coding of information into one of 
M possible signals and transmitting it as a narrow pulse of optical energy 
placed in one of M adjacent time intervals that constitute a word. The nature 
of the technique permits the optical source to be operated at a low, message- 
independent duty-cycle, to extend the lifetime of the device, and furthermore, 
it offers a high degree of immunity to the various noise sources that dominate 
in an optical receiver. However, its high bandwidth requirements make it more 
suited to free-space optical communication links rather than fiber-optic links. 

It has been shown that in the absence of timing errors, the optimal 
decoding scheme for a PPM signalling format is a maximum-likelihood (ML) 
scheme. In this scheme, the receiver compares the number of photocounts in 
each time slot within a word and chooses the slot with the maximum number of 
counts as the slot in which the optical pulse was transmitted. The nature of 
this decoding technique, essentially an energy detection scheme, dictates that 
timing synchronization is of the utmost importance. Timing offsets during 
decoding will cause some of the counts in one slot to be counted as part of the 
photocounts of an adjacent slot. The energy spillover decreases the effective 
signal power, increases the effective noise power and increases the overall 
error probability. 
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The effect of timing errors in optical digital communication systems has 
attracted some attention in recent years [3]. Gagliardi and Haney [4], as well 
as Snyder and Forrester [5], have investigated the use of phase-locked loops 
(PLLs) and the characteristics of the PLL tracking error in optical systems. 
Mengali and Pezzani [6] analyzed the phase error variance of a PLL employed in 
an optical Pulse Amplitude Modulation (PAM) system. 

Lafaw and Gardner [7] have investigated a timing subsystem for optical PPM 
communication links by simulating a receiver employing direct photodetection, 
preprocessing of the detected signal and a PLL for timing synchronization. The 
variance of the relative phase error was studied under varying signal strength 
conditions. 

Further investigations into the system described in [7] have been conduc- 
ted by Chen and Gardner [8], and their results show that the synchronization 
performance improved with increasing signal power and decreasing loop band- 
width. The BER (bit error rate) of the system was also analyzed and compared 
to that for the perfectly synchronized system. In this report, we describe a 
practical optical PPM system that utilizes a semiconductor laser diode at the 
transmitter, an avalanche photodiode (APD) as the photodetecting element in the 
receiver, a signal preprocessor to extract timing information from the detector 
signal and a PLL to achieve timing synchronization. The statistics of the 
phase error at the PLL are analyzed with respect to various system parameters, 
and the results are compared to the computer simulation results reported in [7] 
and [8]. 

In Chapter 2, the overall optical link requirements are discussed and the 
link equations are developed. In Chapter 3, the signalling format and 
detection strategy are described. In chapter 4 the operation of a PLL is 
described and expressions describing the PLL timing performance in terms of the 
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important system parameters are listed. In chapter 5, the various components 
of the experimental set-up in the laboratory are described and our experimental 
results and interpretations are reported in Chapter 6. 
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CHAPTER 2. COMMUNICATION SYSTEM LINK 
2.1 Development of Link Equation 

A free-space optical communication link has similar system components to a 
RF communications link, namely, a transmitter, transmitting and receiving 
antennae and a receiver. The components, however, differ in form. An optical 
transmitter usually consists of a modulator and a light source. The light 
source can be coherent (laser or laser diode) or incoherent (LED). The 
receiver consists of a telescope and a photodetector (PIN photodiode, APD or 
PMT) followed by an amplifier. Figure 2.1 shows the block diagram of such a 
link. 

Consider first the optical transmitter. For spatially coherent laser 
sources, the transmitting telescope is usually a beam expander that is designed 
to reduce the beam divergence. The beam expander may consist of reflecting or 
refracting optics depending on cost, weight and wavelength of the source. 
Because of the high cost of low-loss lens materials, reflecting optics are 
utilized almost exclusively in the middle and far infrared wavelengths. Both 
types are used for near infrared and visible systems. Figure 2.2 shows two 
common beam expander configurations. 

The fundamental lasing mode for a radially symmetric cavity has a Gaussian 
cross section [9] 

I L (x,y,z) = [2P L /itv 2 (z)]exp[(-2(x 2 + y 2 )/v 2 (z)] (1) 

2 

where 1^ = beam intensity (V/m ) 

P L - total power in the beam (V) 

w(z) = half beamwidth at the e intensity point (m). 

For a collimated beam, 
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Figure 2.2. Typical beam expander configurations. 
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w(z) = w o [l + (Xz/Jtv o 2 ) 2 ] 1/2 (2) 

where 

v q « beamvidth at output of beam expander (m) 

X = transmitted wavelength (m) 

z = distance between transmitting aperture and receiver (m). 

In the far-field region, the on-axis signal intensity due to the 
collimated beam is thus given by 

I s (z) = {2w 2 P L T a /X 2 z 2 }/[l + (iw o 2 /Xz) 2 ] (3) 

where 

T a = atmospheric transmittance. 

For the laser diode in our system however, the lasing cavity is not 
symmetric. As a result, the output laser beam is elliptical in shape, as can 
be observed in Figure 2.3. In addition, the output beam is also highly 
divergent, due to the laser diode cavity's extremely small aperture. To 
correct for these geometric aberrations, a laser diode collimating lens is 
employed at the laser diode output so as to obtain a collimated elliptical 
beam, which is then passed through an anamorphic prism pair in order to obtain 
a circular beam. The prism pair expands the minor axis more than the major 
axis, thus yielding a radially symmetric output beam. The circular beam is 
then sent through a beam expander whose output is directed onto the plane of 
the receiving antenna. Figure 2.4 shows the optical transmitter when a laser 
diode is used as the light source. 

Since the Gaussian beam emanating from the coherent light source is 
infinite in extent, Eq.(3) holds for a beam expander of infinite dimensions. 

In practice, the beam expander output is a truncated Gaussian beam. The 
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Gaussian beam approximation is good if the output beam diameter is much smaller 
than the beam expander diameter. Unfortunately, by implementing such a condi- 
tion, a large portion of the beam expander remains unused. This redundancy 
increases the system cost considerably. Klein and Degnan [10] have developed 
the on-axis far-field intensity pattern of a truncated Gaussian beam as 


1 ( 2 ) = 


2hw 2 T P t 
o a Li 


[exp(-( 


w 


-) 2 > 


exp( 


o 


(A) 


where a = radius of transmitter aperture 

b = radius of central obscuration (for Cassegrainian telescope). 

Buck [11] has demonstrated that the maximum intensity in the center of a laser 
beam is obtained by using the largest transmitting aperture available (within 
practical limits) and by expanding the beam so that the intensity drops to 8 
percent at the edge of the output aperture of the beam expander. 

At the optical receiver, the configuration shown in Figure 2.5 is normally 
used as the receiver antenna. The incident beam is focused via an objective 
lens through a field-stop iris onto a collimating lens. It is then passed 
through an interference filter before being focused onto the photodetector. 

The field-stop iris reduces the aperture through which background radiation can 
enter the receiver. The interference filter, which is centered around the 
transmitted wavelength and has a passband of tens of angstroms, further reduces 
the noise level. As such, the received power is given by 


P r (z) = 


n R n T A R 2lla> o T a P L a 2 ,b .2 ,2 

— [exp(-(— ) ) - exp(-(— ) )1 

o o 


x 2 2 

A z 


■ < X/4,c > T a \ \ C R g t P l 


(5) 




Objective 
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Figure 2.5. Simple direct detection optical receiver 
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where Or = receiver efficiency due to receiver optics 

yip = transmitter efficiency due to transmitter optics 
2 

Gr = 4 n Ar/X = equivalent receiver antenna gain 
2 2 2 

Gy = 8n w q /A = equivalent transmitter antenna gain 
Ar = area of receiver objective lens 

Equation (5) is referred to as the link equation, since it gives an expression 
for the received signal power at the photodetector in terms of the 
communication link parameters. 

2.2 Signal and Noise Analysis at the Detector 

The photodiode is a semiconductor device which converts incident photons 
into electrical impulses. For each photon packet that strikes the surface of 
the photodiode, a primary electron-hole pair is created in the bulk region of 
the device. Device parameters can be manipulated during fabrication so as to 
optimize their responsivity to specific ranges of wavelengths of the incident 
radiation. The nature of the operation of photodiodes is well characterized in 
the literature and has been found to obey quantum statistical laws. 

For a PIN photodiode, which has no internal gain, the expected primary 
photon current generated by a single signal photon impinging upon its surface 
is given by 

<i s > = ^ P s /hv = RP S < 6 > 

where 

0 = quantum efficiency of the photodiode 

- 19 

q = electronic charge = 1.601x10 C 
P s = incident signal power 

h = Planck's constant = 6.626x10"^^ W-se c/Hz 
v = frequency of the incident radiation 
R = responsivity 
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An avalanche photodiode (APD), has an internal gain mechanism that can be 
characterized as a random process- For an APD, the expected primary photon 
current due to a single photon is given by 

<i s > - mP s <G>/hv (7) 

where <G> = average gain. 

There are several noise sources at the receiver of an optical communica- 
tion link. To optimize the receiver design, the various noise sources have to 
be investigated and receiver components carefully selected so as to maintain a 
given signal to noise ratio (SNR). The power SNR at the output of an optical 
receiver is defined as 

__ total signal power from signal photocurren t 

photodiode noise power plus thermal noise power 

The photodiode noise consists of incident background noise as well as shot 
noise, which arises from the statistical nature of the current generation 
process in the photodiode. An APD exhibits avalanche multiplication noise as 
well, which is dependent on the mean-squared gain. The overall noise model for 
the photodetector is depicted in Figure 2.6. It can be shown that [12]: 

2 

<Gi > 

s 

2q(i s +i b +i d )G 2 F(G)B + 2qi 1 B + 4kgTB/R L 


SNR = 


( 8 ) 




i - Signal current 

i ds - Gain-independent dark current 

i t - Preamp equivalent thermal ncise current 

i b - Current due to background radiation 

i - Gain-dependent dark current 

G - Average detector gain 

R - Amplifier input resistance 

C - Detector capacitance 

d 

C - Amplifier capacitance 

3 


Figure 2.6. Direct detection receiver circuit model. 
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where 

i = current due to signal power (A) 

= current due to background noise power, P^ (A) 
ij = dark current of APD (A) 
i^ = surface leakage current of APD (A) 

G = <G> = average gain 

B = detector bandwidth 

k D = Boltzmann's constant = 1.379xl0~ 33 w/K-Hz 

D 

T = absolute temperature (°K) 

= APD load resistance (2) 

F(G) = excess noise factor of APD. 

The excess noise factor is given by the relation 

F(G) = k eff G + (l-k eff )(2-l/G) 

where k ,, = effective ionization ratio of holes and electrons in the material, 
eff 

For Si, k ^ = 0.02, and the following approximations can be made: 

F(G) = G°* 3 10 < G < 150 

F(G) = 0.09G 0 ’ 9 150 < G < 1000 

Photodiodes can be designed to operate in either the shot noise or thermal 
noise limited regions. Thermal noise limited detection occurs when the 
internal gain of the detector and the optical signal power are low. In this 
case, thermal noise dominates the detection process. Shot noise limited 
detection occurs when the internal gain of the detector is high so that the 
shot noise dominates. The SNR expression can be reduced to simpler forms for 
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CHAPTER 3. SIGNALING FORMAT AND ERROR RATE ANALYSIS 
3.1 Pulse Position Modulation 

The detection strategy employed in the receiver is dictated by the 
signaling technique employed. Once this has been decided upon, the performance 
of the system is evaluated in terms of the system's error rate. The error 
rate can be expressed in terms of the SNR, vhich was derived in the previous 
chapter. 

The signaling format under consideration is M-ary pulse position 
modulation (PPM). In PPM, a single pulse is sent in one of M time slots of 
duration T g seconds. Figure 3.1 depicts the relevant parameters of a general 
PPM scheme. The PPM scheme is a highly efficient form of block encoding, since 
log 2 M bits of information can be transmitted via a single pulse of light. Due 
to the large difference in energy levels between the signal interval and the 
M-l noise intervals within a word frame, a particularly simple detection 
strategy can be employed at the receiver for a PPM system. Figure 3.2 shows a 
block diagram of a PPM block decoder. The decoder integrates the photodetector 
output current over each of the M time intervals and picks the largest of the M 
values as the slot in which a signal pulse was sent. In essence, this strategy 
is an energy-detection scheme. 


3.2 Error Rate Analysis 

The performance of block encoded systems is indicated by the average 
probability of making a word error, PWE. This is given by [ 13 ] 


PWE 


1 


exp[-(K s+ MK b )] 

M 


- lPos(k,K s+ K b ) 
k=l 


k-1 

x [ } Pos(j ,K b ) ] M_1 (j^) [(l+a) M -l] 

3=0 


(ID 



Number of time slots per word = 
Number of bits per word 
Total time per word 

bits per word 

Data Rate = 

time per word 

Duty Cycle 


M 


log 2 M 

MT 


log 2 M 
MT 


p 

MT S 



Figure 3.1. Pulse Position Modulation (PPM) scheme 










where 



k! j K J/t! 
t=0 


Pos(k,y) exp(-y) 

K = average signal count per word 
s 

= average background count per slot. 

The probability of bit error, PBE is related to the PWE by 

PBE = j (j^-) PWE 


( 12 ) 


Note that for thermal noise limited detection, with M = 2, [9] 

PBE = i [1-erf (/SNR/2)] 

where SNR is given by Eq. (9). The above relations can be written in terms of 
the parameters listed previously. Assuming that rectangular pulses are 
transmitted, then 


where A g is the expected photocount rate and is given by 
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where f| = photodetector quantum efficiency 
P g = incident signal power (W) 
h = Planck's constant = 6.62xlO - ^ tf-sec/Hz 
\> = optical frequency (Hz). 

Similarly, 

h ■ Vs 

= tiP^/hv (sec - *) 

where 

P^ = incident background noise power (W). 

Figures 3.3 and 3.4 show how the PWE for an M-ary PPM system varies with 
K g , and M. These are some of the critical parameters which affect the 
performance of the PPM communication link. 

Another key factor which has a bearing on the performance of a PPM link is 
that of timing accuracy. Timing offsets during the decoding process cause 
energy spillover from a time slot into its adjacent time slot, resulting in 
intersymbol interference. The spillover is tantamount to a decrease in the 
effective signal power, increased effective background noise and increased 
probability of error. Figure 3.5 shows how the PWE varies with timing error. 


PWE 



Figure 3.3. Comparison of the exact error rate and the 
union bound estimate for a 4-ary PPM system 
(From [8], p.28) 
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CHAPTER 4. THEORETICAL ANALYSIS OF THE PHASE-LOCKED LOOP (PLL) 

In Chapter 3, we discussed the stringent timing requirements of an optical 
PPM communications link, and observed how the PWE is affected by a fixed timing 
error. For receivers using dynamic phase synchronization circuits such as a 
PLL, the receiver timing error is random, in general. The statistics of the 
timing error can be expressed in terms of the statistics of the phase error of 
the PLL employed in the receiver. PLL parameters, which have been extensively 
studied [14], can thus be used to indirectly determine the timing error 
statistics . 

4.1 Mathematical Model of PLL 

Figure 4.1 shows the block diagram of a PLL. Figure 4.2 is the equivalent 
baseband phase model of the PLL. A phase comparator determines the phase 
difference between its two inputs and generates an error voltage which is 
proportional to this difference. The error voltage is filtered and fed into a 
voltage-controlled oscillator (VCO), which adjusts its frequency to minimize 
the phase difference between the VCO output and the incoming signal. The PLL 
achieves "lock" when the VCO frequency is equal to the incoming signal 
frequency. 

For the PLL analysis, we assume that the input to the PLL consists of a 
sinusoidal signal, V^t), as well as an additive noise term, n(t). V^t) is 
assumed to have amplitude, A^, phase, and radian frequency, w. . The phase 
detector of the PLL is viewed as a perfect multiplier with gain, K m , and the 
VCO is assumed to have a gain constant, K- v , and a free-running frequency, (*) q . 
Hence, the frequency detuning, G = | - oo q | . From the baseband model of the 

PLL, it can be shown that the phase error between the incoming signal and the 
VCO output signal, 9 , satisfies the following nonlinear stochastic 
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differential equation [15]: 

d9 (t) A 

-jpf = 2 - K y f ( t ) * sin © e ( t ) + n(t) cos(w.t + © Q ( t ) ) ] (13) 

where f(t) is the impulse response of the low-pass filter. Under locked 
conditions, 0^(0 is small, so sin[© e (t)] ~ 9 e (t). We can thus approximate the 
above model with a linear PLL model, depicted in Figure 4.3. The transfer 
function of the linear PLL model is 


H(s) 


e 0 ( s) 

e.(s) 


KF(S) 

s+KF(s) 


(14) 


where H(s) is the Laplace transform of the PLL system impulse response, F(s) is 
the Laplace transform of the loop low-pass filter, and K = K^K^/2. 

Note that H(s) is the system transfer function when no noise is present at 
the input. If we choose a first-order loop filter to avoid complexities, then 
the loop transfer function can be written as 

F(s) = a/(s+a) (15) 

where a = 1/RC 

and R and C are the loop filter components. It is now straightforward to show 
that, when the input phase is written as 

9.(t) = w.t + e lt 

where 0^ is the initial input phase, then the Laplace transform of the phase 
error can be written as 

(s+a)(9 + s9^) 

®e^ S ^ - s[s(s+a) + Ka] 
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Applying the final value theorem to Eq.(16), we observe that the steady-state 

phase error, when we let t -» is 

0 = S2/K . 

00 

Thus, it is evident that the steady-state phase error is minimized by 
minimizing the frequency detuning. 

With Eq.(15) as the choice for F(s), the PLL transfer function becomes 

H(s) = (17) 

s +as+Ka 

Some of the other important PLL parameters are oo^, the natural frequency of the 
system, £, the damping ratio, and the hold-in and lock-in ranges. The hold-in 
range, Aoc^, is the range of frequencies over which the loop will remain in lock 
and is given by 

A(a^ = 2KF(0). 

Note that A is the dc loop gain of the system. The capture range, Aoo^, is 
the frequency range over which the loop acquires phase-lock without slipping 
cycles. It has been shown that for a pseudorandom input signal, the capture 
range is approximately equal to the loop noise bandwidth [15]. The 2-sided PLL 
noise bandwidth, (rad/sec), can be defined as 

00 

2 

B T = J | H ( go) | doo for Fourier Transforms 

_ 00 

1 1 

= J* H(s)H(-s)ds for Laplace Transforms. (IB) 

-i°° 

It has been shown that when B^ is small, the phase error of the PLL, 
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© e (t), is a wide-sense stationary (VSS) process (16]. More specifically, 9 g (t) 
is a WSS process if the loop transfer function satisfies the following 
condition: 


H(oo) = 0 for |o>| > n/T 


(19) 


where T is the word period of the M-ary PPM modulation [ 8 j . 
w 

For the above PLL model, using Eqs. (14), (17) and (18), and Refs. [15] 

and [17], we obtain the following parameters: 

. 1/2 


to = (Ka) 
n v 


5 = (a/4K) 

Aoc^ = 2K 


Aw = B. 
c L 


1/2 


PLL natural frequency 
Damping ratio 
Hold-in range 
Capture range 


where 


B L = 


K. , a .1/2 
8 4K-a 


K|/x - i/y ] 
4/z 


if a < 4K 


if a > 4K 


( 20 ) 


and x = [2aK - a + z]/2 

y 

z 
K 

a = 1/RC. 


[ 2aK - a - z]/2 
1/2 

= a(a(a - 4K) J 1 ^ 


KK/2 
v m 


4.2 Signal Preprocessing 


It was shown by Chen and Gardner |8] that when the PPM pulses are 
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rectangular, the photodetector output signal does not have a frequency 
component at the slot frequency. Since PLLs require a strong frequency 
component to achieve phase-lock, the photodetector output must be suitably 
processed to extract frequency information. Since any linear operation on the 
signal would only have the effect of multiplying the frequency spectrum of the 
signal by the magnitude square of the transfer function of the linear operator, 
no frequency component is introduced at the slot frequency using a linear 
preprocessor. Consequently, nonlinear processing must be used. One technique 
is to pass the signal through a linear filter followed by a squaring loop [18]. 
We have selected a differentiator, followed by a squarer. The differentiator 
generates alternate positive and negative going impulses at each level 
transition. The squarer ensures that only positive impulses appear at the slot 
frequency. A block diagram of the signal preprocessing circuit is shown in 
Figure 4.4. 


4.3 PLL Performance with Signal Preprocessing 

Assuming that the phase error 0^(0 of the PLL is approximately Gaussian 

distributed, the PLL synchronized system performance is best indicated by the 

2 2 
phase error variance, c r^. Derivation of the exact expression for <x^ is 

very involved, but Chen and Gardner [8] have developed the following 

approximation for the phase error variance: 


M_ A, 

" K ( o> > Y 
s s 


( 21 ) 


2T JJ P(oo' )P*(oo" )P*(w' -oo")G(oo' )G(oo-oo' )G* (oo" )G*( oo-oo" )doo' doo M 

S _ ao 

I P(w)G( oo)*P( w)G(«) | 2 


00=00 


_ ao 
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where K 

s 

M 


Y 

P(«) 

G(w) 


signal count per word 
PPM order 

PLL closed loop bandwidth (radians/sec) 
slot frequency (radians/sec) 
dimensionless parameter which depends on 
Fourier transform of received signal pulse shape 
D(oo)*K(oo) 


D(oo) = Fourier transform of impulse response of APD 

K(oo) = transfer function of linear filter in preprocessor. 

Equation (21) indicates that the phase error variance is inversely 

proportional to the received signal power and directly proportional to the PLL 

bandwidth. Figures 4.5 and 4.6 show how the phase error variance is affected 

by varying the signal count, and the PPM order, M, respectively. Recall 

from Chapter 3, that the critical parameter being investigated here is the 

2 

timing error, e = A/T . The timing error variance, cr for a PLL synchronized 
M-ary PPM system, can be related to the PLL phase error variance by 


2 2 2 
= < /(2rt) 


( 22 ) 


The above relation can be applied to the error rate analysis of Chapter 3 to 
show how the timing error affects the overall error rate of the PPM optical 
communication link. Figure 4.7 shows that as the loop bandwidth increases, the 
PBE of the system increases. 



PHASE ERROR VARIANCE (rad 



Figure 4.5. Computer simulated phase error variance versus the signal count 

rate for a 4-ary PPM receiver. The solid curve is the theoretical 
variance predicted by Equation (21). (From [8], p. 24) 



PHASE ERROR VARIANCE (rod 2 ) 


Figure 4. 



. Computer simulated phase error variance versus the order of PPM 
for an M-ary PPM receiver. Again, the solid line is the theo- 
retical variance predicted by Equation (21). (From [8], p. 25) 



PBE 


Figure 4. 



. Probability of bit error (PBE) of a PLL synchronized 4-ary PPM 
system versus signal count/word for various loop bandwidths. 

Y = 10, K, = 1. (From [8], p. 34) 
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CHAPTER 5. DESCRIPTION OP EXPERIMENTAL SYSTEM 

The system being considered here is an attempt to simulate in the 
laboratory a satellite-based optical communication link with very compact 
transmitter and receiver components. We have utilized a semiconductor laser 
diode as the light source, and an avalanche photodiode (APD) as the 
photodetecting element in the receiver. A pseudorandom code is converted into 
the PPM format and treated as the input signal to the laser diode. The 
pseudo-random code being used is a 16-stage maximal length sequence. The APD 
is followed by a transimpedance amplifier, which offers the best dynamic range 
and sensitivity for a given amplifier noise level among the various receiver 
front-end designs [19]. The preprocessing circuit, as was mentioned in Chapter 
4, consists of a differentiator and a squarer. The subsequent signal is then 
fed into the PLL. In this chapter, we synthesize the information governing the 
operation of the various system subsections described in the previous chapters 
with the aim of selecting the overall system parameters. 

5.1 Transmitter Design 

The transmitter consists of three main components. The front end of the 
transmitter is the maximal length code sequence generator. For the given set- 
up, a 16-stage binary shift register maximal-length sequence is generated, 
which yields a code length of 65535 bits with a variable clock rate of .5, .33 
or .25 MHz. The second part of the transmitter converts the PN code into an 
M-ary PPM signal. Three choices of M are available (4,8,16, corresponding to 
2-, 3- and 4-bit words, respectively). The conversion process first inputs the 
selected bits into an up-down counter, which counts at M times the selected 
clock rate. This is referred to as the "chip" rate. The counter emits a pulse 
when the value in the counter register reaches zero. Hence, only one pulse 
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appears in a single time slot within a word consisting of M time-slots. Table 

5.1 (p.40) depicts the PPM parameters, described in Figure 3.1, for our system. 
Figure 5.1 (p.41) is the schematic of the PPM PN-code generator. This PPM 
signal is then fed as input to the third part of the transmitter, the laser 
diode driver circuit. 

The laser diode driver circuit, whose schematic is seen in Figure 

5.2 (p.42), includes a slow starter section to protect the diode from sudden 
power surges as well as a feedback section. The feedback mechanism maintains 
the laser diode at a fixed bias current just below the threshold current, and 
also maintains a constant output power level, as the laser diode deteriorates 
with time [20]. Figure 5.3 (p.43) shows the actual laser diode driver circuit. 
In this circuit, a SHARP LT024 MD laser diode, capable of a maximum optical 
output power of 30 mW at a wavelength of 780 nm, is employed in conjunction 
with a laser diode collimating lens supplied by D.O. Industries. The laser 
diode operates at an average CW output power which will maintain specific 
signal count levels at all PPM orders. The laser diode is biased just below 
the threshold for a digital input 'O', so as to enable rapid turnoff whenever 
the input signal makes a high-low transition. The laser diode parameters and 
the driver circuit operating point parameters, which are computed in the 
appendix, are listed in Table 5.2 (p.44). The laser diode and the diode 
collimating lens are mounted in a temperature stabilized enclosure, which is 
shown in Fig. 5.4 (p.43). Figure 5.5 (p.45) depicts the schematic of the 
temperature control circuit [21]. Due to the close proximity of transmitter 
and receiver in the laboratory, we have omitted using any other optical 
components as part of the transmitter, which would normally be required to 
reduce the beam divergence. Instead, the beam is passed through a cascade of 
several neutral density (ND) filters to simulate atmospheric attenuation 














Figure 5.1. Schematic of PN code PPM signal generator. To start the code, 
the switch is toggled from position A to B, then back to A. 















Figure 5.2. Laser diode driver circuit schematic. The circuit is activated 
by flipping the switch from A to B. 















Table 5.2 


Laser diode and driver circuit operating point parameters 


Laser diode specifications 


e // - 

CO 

0 

= FWHM angle of laser beam measured parallel to junction plane 

0a. = 

31.2 ° 

= FWHM angle of laser beam measured perpendicular to junction plane 

' th = 

43.1 mA 

= threshold current 

T1 = 

0.74 W/A 

AP 

= differential efficiency - 1 ^ | 379 ' 

X = 

785 nm 

= wavelength 


Laser diode driver circuit operating point parameters 


M 

P ave< mW > 

4 

2.65 

8 

2.12 

16 

1.71 


1 bias = 39 5 mA 

P 


0.44 mW 



T| - Fenwal FB3IJI-WC Thermistor 

TEC - Melcor CP 1.4-31 -IOL Thermoelectric Cooler 
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Figure b.b. Schematic of thermoelec t r ically-cooled temperature control 
circuit for the laser diode* 
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and space loss. The filters can be varied to effect different received signal 
strengths. In our system, a tungsten lamp, with a variable polarizing filter 
placed in front of it, is treated as a white light source for incident 
background radiation (noise). 

5.2 Receiver Design 

At the receiver, a narrow-band interference filter, centered at 780 nm, 
with a half-maximum bandwidth of 10.8 nm and a peak transmission of 49.5%, is 
used. Both the signal beam and the noise beam are incident upon the filter. 

The composite filtered beam is then focused onto the detector surface. The 
receiver fixture is designed so that the APD surface and a power meter can be 
alternatively positioned at the focused spot, which enables the incident light 
power to be measured. This arrangement is shown in Figure 5.6. 

Figure 5.7 depicts the schematic of the receiver. The front end of the 
receiver consists of a Fujitsu FPD08R32WS Si APD, connected to a National 
Semiconductor LH0032CG wideband op amp, which is configured as a transimpedence 
amplifier. Both devices are contained within a thermally insulated box, the 
inside of which is kept at a constant temperature of 26°C by means of a 
temperature control circuit very similar to the one employed in the 
transmitter. The reason for the temperature control is two-fold. The fact 
that the APD gain varies significantly with temperature fluctuations makes 
temperature stabilization mandatory. In addition, the wideband op amp 
generates a considerable amount of heat, thus making thermal noise the dominant 
noise source at the receiver. 

The gain control amplifier is located at the output of the transimpedance 
amplifier. The receiver signal can thus be suitably amplified to ensure that 
the input signal to the PLL has a constant amplitude whenever the PPM order is 
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Figure 5.6. 


Flipping mirror arrangement at receiver facilitates 
measurement of optical power incident on photodetector. 


A 




Figure 5.7. Schematic of optical receiver. 
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changed. Changing the PPM order causes the laser diode extinction ratio to 
vary. Hence, the gain control is required so as to make an unbiased comparison 
of the effects of varying the PPM order on the performance of the receiver. 

The signal preprocessing segment of the receiver consists of a simple R-C 
differentiator followed by a squaring device, a Motorola MC 1595L wideband 
multiplier. The squaring device has a gain constant which is determined by the 
variable resistors Rl, R2 and R3 in Figure 5.7. 

The timing extraction section of the receiver consists of a Signetics 
NE 564 PLL as configured in Figure 5.7. The PLL generates an error voltage 
proportional to the phase error, 

V') - W') - V [)l 

where K is the phase detector gain and 0. and 0 are the phases of the input 
m 10 

signal and VCO signal, respectively. On the NE 564, is measured between 
pins 4 and 5 [22], The differential signal between these two pins is buffered 
and then translated into a single-ended signal to facilitate data acquisition. 
The complete receiver board is shown in Figure 5.8. 

The loop constants for various loop bandwidths and VCO free-running 
frequencies were determined experimentally and are tabulated in Table 5.3. The 
loop bandwidths, B^, were computed using Eqn. (20). From the last column in 
Table 5.3, it is obvious that Eqn. (19) is satisfied, and hence, the phase 
error is deemed to be a VSS process. 
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Figure 5.8. Receiver circuitry, comprising gain-control amplifier, 
differentiator, squarer, PLL, buffer amplifiers and 
temperature-control circuit. 





6800 2.348 0.3318 7.261 23.173 0.0087 

2200 3.941 0.2713 12.217 48.882 0.0122 

4.000 4700 3.940 0.2826 12.217 33.032 0.0083 

6800 3.938 0.2606 11.719 26.027 0.0065 


input frequency 

loop bandwidth capacitor 
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5.3 Data Acquisition 

Figure 5.9 is a block diagram of the experimental setup. The actual 
laboratory setup is seen in Figure 5.10. The PLL phase error signal is measured 
using a digital voltmeter. The experimental results are based upon a total of 
135 sets of data, with each data set comprising 15 data samples. The acquired 
data is processed and presented in the following chapter. 



ND 



Figure 5.9. Block diagram of experimental setup. 




V r ' jy A}/ " 'y. 


ri 

mental setup of 

th 

e optical communications 1 

in 

tr 

ansmitter is at 

r i 

ght foreground, the white 

no 

ce 

is at center ba 

ck 

ground, and the receiver 1 

ie 

0 

f center, in the 

f 

oreground . 






55 


CHAPTER 6. EXPERIMENTAL RESULTS 

In Chapter 4, it was mentioned that for our system to function properly, a 
slot- frequency component of the detected signal had to be present for the PLL 
to lock onto. Figure 6.1 shows the effects of the receiver preprocessing 
section on the detected signal. 

The power spectrum of the detected PPM signal shows no slot frequency 
component at 2 MHz (Figure 6.1(b)). A small 2 MHz signal component is seen at 
the differentiator output (Figure 6.1(c)) and the nonlinear operation 
(squaring) increases the signal strength of the 2 MHz component by about 15 dB 
(Figure 6.1(d)). Hence, the signal preprocessing circuit achieves its 
objective of evoking a frequency component for the PLL to acquire lock. 

2 

Figures 6.2, 6.3, 6.4 and 6.5 show how the PLL timing error variance, a , 

changes with different system parameters. Figure 6.2(a) presents the variance 

versus signal power plots for differing background noise levels. Figure 6.2(b) 

is the linear least-squares (LLS) fit to the data presented in Figure 6.1(a). 

The plots indicate that the timing error variance increases with signal power. 

2 

The dependence of a on the background noise is inconsistent, though, 

increasing when P n is increased from 0 to 75 nW, but decreasing to the original 
o 

level when P D is further increased to 134 nW. 
o 

2 

Figure 6.3(a) plots <j versus the signal power for different PPM orders, 

V 

and Fig. 6. 3(b) is the LLS fit to the data in Fig. 6. 3(a). Again, it is observed 
that the timing error variance increases with signal power. Contrary to the 
behavior predicted by Eq.(21), it appears that the variance decreases when the 
PPM order is increased. However, it should be noted from Fig. 6.3(a) that, for 
the M=4 and M=S cases, some widely varying data points occur, which may be the 
cause of the disagreement between theoretical prediction and the presented 


data. 
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Figure 6.1. Effects of signal preprocessing on detected signal. In (a), the 
oscilloscope screen shows the PPM signal for the M=4 case, (b) 
is the frequency power spectrum of the signal in (a). The slot 
frequency is 2 MHz, but no frequency component at 2 MHz is seen 
in (b). (c) shows the frequency spectrum of the differentiated 

signal. A frequency component is seen at 2 MHz. In (d), the 
frequency spectrum of the signal at the squarer output is seen 
to have a stronger 2 MHz component than in (c). 
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2 

Figure 6.4(a) shows a versus signal power for different PLL loop 

bandwidths, B^. Some inconsistencies are noted in Fig. 6.4(b), the LLS fit to 

the above data. For the two smaller settings, the variance is observed to 

decrease with increasing signal power. The converse is true with the largest 

B^ setting. A further inconsistency is the observation that, when B^ is 

2 

increased from its smallest value, a decreases. However, when B. is further 

£ r* 

2 

increased, <x increases. Figure 6.5 shows how the timing error variance 

changes with PPM order, at fixed signal and background noise levels. The plot 

2 

shows a marked increase in cr ^ when M is increased. 

£ 

The experimental results obtained thus far have not shown a high degree of 

correlation with Eq.(21). Table 6.1 lists the gradients and y-axis intercepts 

of each of the 9 LLS estimate curves plotted in Figures. 6.2, 6.3 and 6.4. In 

2 

an attempt to estimate the relation between the timing error variance, a and 
the signal power, P g , the following model is assumed: 


2 

a 


£ 


P 

s 


S 


Eq.(21) can be written equivalently as 


2 M , B L, 

% = r Y 

s s 


2 2 

= (2 It; a\ 


where K 


hP T \<G> 
s s 

he 


2 

u 


he My 
2 



(23) 


s 


hX<G>(2n) 


(24) 
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Table 6.1 Gradients and Y-axis intercepts of the 9 linear least- 
squares fit curves. 


I 

I 

I 


Figure 

Curve 

log a 

(Y-axis intercept) 

Y 

(slope) 

a 

6.2 

P s = 0 nW 
P s = 65 nW 
P s = 1 34 nW 

- 3.106355 

- 1.925988 

- 3.185269 

1.196060 
1.278143 
1. 157769 

0.0448 
0. 1 457 
0.0414 

6.3 

M = 4 
M = 8 
M * 16 

4.629758 
0.3630204 
- 17.38993 

1.764760 

1.410833 

0.163823 

1 02.489 
1.4377 
2.803E-9 

6.4 

B L =2 1 0 1 5 Hz 
B L = 24576 Hz 
B l = 37975 Hz 

-22.44087 

-34.51324 

-12.91279 

-0.084518 

-0.974457 

0.5396607 

I.795E-I0 

I.0259E-15 

2.4663E-6 
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hcB^My 

where a = *— 

nA<G>(2Ttr 

From Eq. (23) 

2 

log ( (J^) = log a + S log P g 

Using the entries in Table 6.1, we get 

2 

log ( a^) = <log a> + <S> log P g 
s (Y-intercept ) + (slope) log 

where <log a> and <8> are the LLS estimate of log a and S respectively. 

The objective of this exercise is to estimate y> the dimensionless constant 
in Eq. (21), which is a function of the pulse shape and the filter impulse 
response. For the system being considered here, y has to also include the 
effects of the gain amplifier and the squaring device. 


The linear least-square estimates of y> are given in Table 6.2. 


Figure 

Curve 

6.2 

P s = 0 nW 
P s = 65 nW 
P s = 1 34 nW 

6.3 

M = 4 
M = 8 
M - 16 

6.4 

B l =2 1 0 1 5 Hz 
B l = 24576 Hz 
B l = 37975 Hz 
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CHAPTER 7. CONCLUSIONS 


The results obtained from the experiment described in this paper were 
inconclusive, as far as verifying the theoretical model is concerned. Several 
factors could have contributed to the supposedly erroneous experimental 
results. The biggest factor affecting the results is probably the data 
acquisition technique. A computer-controlled data acquisition system with a 
sampling rate of at least 100,000 samples per second is required for the data 
samples to be stochastically independent. Hand-held measurements using a 
digital voltmeter, such as were performed during this experiment, are highly 
inappropriate for this purpose and introduce large errors in the data. The 
experimental observations indicate that for the most part, the timing error 
variance increases with signal strength. This would imply that the operating 
point of the system lies in a region where noise, rather than the signal, 
dominates. One possible source of the noise may be electromagnetic 
interference (EMI) due to lack of adequate shielding. 

Further refinements can be made to the above described experiment. 

Firstly, a computerized data acquisition procedure must be adopted. 

Furthermore, adequate shielding must be provided to the system so as to 
minimize the effect of EMI. In addition to these measures, it is also 
suggested that a more sensitive APD be used, so that the system can be operated 
at received signal strengths closer to the photon counting region. Once the 
theoretical predictions have been confirmed, the experiment may be taken a step 
further by constructing a PPM decoder, so that the system error rate can be 
investigated . 
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SELECTION OF OPERATING POINT PARAMETERS FOR THE LASER DIODE 

The PPM transmitter has the following parameters: 


M 

T s (nsec) 

4 

500 

8 

375 

16 

250 


When investigating the effects of varying PPM order on phase error 

2 

variance, we have to fix the signal count, Ks, and compute at different 
values of M. From Chen and Gardner [8], we have 


2 M B L 

ff <|> = K~ 

s s 


K 


s 



T 

s 


where P s = average incident power 

If we normalize K g , (ie, = 1), we have the following: 
M = 4; 

K s = ( F^ )(500 nsec) P s4 = 1 
= 125 nsec > 4 P s4 = 


1 




M = 8 
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K 


s 


<B )<375 nsec) P s8 


1 


M = 16; 


(1L 125 nsec) 3 P s8 


K = 


(1L)(250 nsec) P sl6 


= 1 


" ^ 125 nsec) 2 P sl6 


= 1 


Hence, 4P g4 = 3P g g = 2P sl6< 

Let Pq = laser diode output power for digital '0' signal 

P^ = laser diode output power for digital '1' signal 

Thus, if we choose P, = P for M = 4, to maintain a fixed K , we require 

J. 5 

that 


M 

P 1 

4 

P 

8 

4/3 P 

16 

2 P 


Let P = 15mW, so that for M=16, doesn't exceed 30 mW, the rated maximum 
power of the LD024MD. This yields 
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M 

P^mV) 

4 

15 

8 

20 

16 

30 


P 

ave 


P 0 (M-1) + Pj 
M 


(i) 


If the extinction is to be at least 50 (arbitrary), then choose Pq 
U sing equation (i), the following operating points are obtained: 


M 

P (mW) 

ave 

4 

3.98 

8 

2.76 

16 

2.16 


.3 mW 


Note that the above are maximum rated values. 
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